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TRANSIENT CRATER MOTIONS: SATURATED SAND CENTRIFUGE EXPERIMENTS 

Schmidt, R. M., Boeing Aerospace Company, M/S 13-20, Seattle, WA 93124 

The kinematics of crater formation are being investigated using a quarter-space test bed 
on a centrifuge. This technique (1) allows high-speed movie coverage of target deforma- 
tions as seen in cross section through a 2-inch Ple/.iglas window along a diametrical 
symmetry boundary. The data for the eight shots performed to date with the recently 
installed Fastax II-100 rotating-prism camera are shown in the attached table. Two 
objectives of these tests were to dynamically record the maximum depth of the transient 
crater and to record its collapse. Preliminary interpretation of the film records for shot 
307-QXE is shown in the attached figure. This shot was conducted in nearly-saturated 
moderately-dense fine-grained Ottawa ’’Banding" sand. 

For this target material, crater floor uplift ("collapse") occurs on the same time scale as 
the radial motion of the ejecta plume. Because of the impairment of the optical path due 
to detonation products plating onto the window, the floor uplifting sequence could not be 
resolved in some of the intermediate frames of the work print. It appears to be complete 
by frame 21 and maybe earlier. At this time the ejecta plume is in the vicinity of the 
final apparent crater radius and the crater profile is in agreement with the final apparent 
crater measured after the shot with the profiiometer. Crater volume displays power-law 
growth until maximum transient depth occurs (frame 4). More interestingly crater radius 
shows power-law growth all the way out to final crater radius (frame 21). Both time 
exponents, interpreted with the assumption of a unique coupling parameter, 6a/(3+<X) for 
volume and 2a/(3+ot) for radius, are in good agreement with a value of the gravity 
exponent a equal to 0.65. This same value was determined for crater growth in water (2) 
(3) and may be applicable for all non-porous geological targets. 

As an expedient when determining the maximum charg' that would not spall or 
puncture the window, the explosive charges were built up Dupont Deta-Sheet. This 
is PETN with approximately 25% inert plastic filler, which is presumed to cause much of 
the smoke and detonation debris products. Pure PETN is now under consideration for 
future shots to improve resolution by reduction of smoke and opaque detonation products. 
However, since the Deta-Sheet detonation products deposit on the window as it is exposed 
to the transient crater, a permanent record of the envelope of maximum crater depth 
versus range is formed. Preliminary comparisons of the detonation product profile with 
the corresponding final apparent crater profile along with viewing the films indicate that 
there was no significant difference in the transient crater phenomena over a variation of 
approximately an order of magnitude in sand grain size. This observation needs further 
examination when the original film records are digitized. 

For seven of the eight shots, the water table was nominally at 1.05 cm below the original 
ground zero. This was to correspond to the approximate scale of the 500-ton PRAIRIE 
FLAT explosive crater in the field, which was modeled in the tangent-above shot 306-QXE 
(see ref. 4). 

(1) Schmidt, R. M. <Sc Piekutowski, A. 3. (1983) Lunar <5c Planetary Science XIV, p. 668. 

(2) Schmidt, R«, M. & Holsappie, K. A. (1982) Estimates of crater size for large-body 

impact: Gravity-scaling results. Geol. Soc. of America SP 190, p. 93. 

(3) Schmidt, R. M. (1983) Trans. Am. Geophys. U. (EOS) Vol. 64, No. 45, p. 747. 

(4) Schmidt, R« M., Fragaszy, R. 3. <Sc Holsappie, K. A. (1981) Proc. of Seventh Int. Symp. 

on Military Applications of Blast Simulation, Medicine Hat, Canada, p. 4.2-1. 
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TRANSIENT CRATER MOTIONS* SATURATED SAND CENTRIFUGE EXPERIMENTS 

Schmidt R. M. 

QUARTER-SPACE DATA TABLE 


EXPLOSIVE MATERIAL CRATER rr-GRCUPS 


Shot No. 

Gravity Mass Q Density Type 
(G) (gm) (C*rg/g) (g/cc) - 

Density 

(g/cc) 

V r h V 
(cm 3 ) (cm) (cm) (cm) 

I* 




y 

306-QXE * 399 

13.0 

■Tfl, $♦ 

1.48 

W8S' 

2.06 

232 

8.84 

2.08 

- 

ItOxlO " 1 

36.8 

4.67 

1.13 

- 

307-QXE 

400 

4.45 

365 

1.48 

WBS 

105 

289 

8.86 

2.77 

4.68 

1.47*10** 

134 

6.85 

2.14 

3.62 

319-QXE“403 

5.49 

415 s 

1.45 

WBS 

2.06 

232 

7.45 

195 

- 

1.45*10** 

87.1 

5J7 

2.13 

- 

324-QXE 

398 

4.23 

42 JS 

1.45 

was 

2.06 

202 

7.21 

2.60 

- 

1.31*10-* 

98.4 

5.67 

2.05 

- 

J26-QXE 

398 

4.14 

42.5 

1.45 

WBS 

2.06 

212 

7.75 

140 

- 

1.30*10“* 

105 

6.14 

j.90 

- 

328-QXE 

401 

4.25 

42 JS 

1.45 

WFS 2 

2.11 

193 

7.41 

2.17 

3.64 

1 .32*1(T* 

95.8 

5.87 

1.72 

2.88 

329-QXE 

100 

4.14 

42.5 

1.45 

WBS 

2.06 

392 

9.76 

325 

6.55 

3.27x10-* 

195 

7.73 

2.58 

5.19 

330-QXE 

100 

4.19 

42.5 

1.45 

WFF 3 

1.97 

320 

8.04 

3.4S 

5.90 

3.29x10-* 

150 

625 

2.71 

4.59 


* Maximum Transient Crater Depth ♦* Tangent Above - ail others Tangent Below «* Water Table * 0.70cm 
' WBS - Wet Banding Sand (212m) 2 WFS • Wet Flintshot Sand (8C0/*) 3 WFF « Wet F-140 Sand (106/*) 

4 Green Deta-Sheet (Type C) 8 Red Deta-Sheet (Type D) 


SHOT 307-QXE TRANSIENT CRATER FORMATION AOOG 
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Cratarfng Flow Fields: A General Form and th« Z Model 
Kevin Houma, Shock Phyilcs and Applied Mechanic*, MS 13-20, Booing Aorotpaeo, Soattlo, WA, 96124. 

The nature of material motlani is often fundamental in discussions of Impact-related processes. 
An analytic 'description of cratering flow fields, known as the Z model, was developed by Maxwell and Seifert 
(1974). Although originally developed for near-surface explosions. It has also been used In calculation 
of Impact ejecta velocities, transient craters and mslt generation (e.g., Schultz and Gault, '1979; Croft, 
1980; Grieve and Cintala, 1981). While it reasonably, describes the flow field In some cases (e.g. Orphal 
at at, 1977), In others it does not. A more general' expression' for flow velocity Is developed here which, 
in certain instances, reduces to' a Z model. 


Z Model 

The Z model was based on a series of numerical simulations of cratering events. Briefly, the radial 
velocity of material behind the outgoing shock wave was observed to vary as a power of the radial position 
r. Discussions of the Z .model often include a time dependent term so that 

r - A(t) r“ z ' (1) 

where A(t) describes the time dependence of the flow. After the passage of the shock, the material 
density was found to be nearly constant. The observed Incompressibility can be used with eg (1) to 
determine the tangential velocity component (e.g. Maxwell, 1977): 

rtf - f (Z-2) *in(tf)/|1+coi(tf)j (2) 

where tf is the angle measured from the downward axis (tf*0 vertically down and tf-tr/2 along the target 
surface). Note that material flows upward toward the surface when Z>2 (tf>0) and downward away from 
the surface when Z<2 (tf<0). Radial flow occurs when Z*2. 

Once A(t) Is specified, egs (1) and (2) determine, among other things, the temporal growth of the 
crater (Maxwell, 1977). A case often adopted in the literature is that of steady flow, 'where A is independent 
of time. In this case, the radius, R c (t) of the expanding crater is 

R c (t) « t 'Az +1 > (3) 

Measurements of crater growth in a variety of materials show that, indeed, crater radius grows as a power 
of time (e.g. Schmidt, 1983,1984; Holsapple, 1984). However, for some materials, steady flow cannot 
hold. For example impacts and explosions in water and wet sand targets (Schmidt, 1984) imply R c grows 
as t 0Ja , corresponding to Z-1.8. But, as noted above, the flow field is rotational away from the surface 
for Z<2. Consequently, a steady flow Z model cannot apply to these cases. 

Austin at c/{ 1980) hove compared *numerical simulations of impact in piasticene clay to predictions 
based on the Z model. They found It necessary to modify the basic form of the Z model. The flow 
center had to be located beneath the surface and both A and Z increased with time. In particular, Z 
was somewhat less than 2 early on end steadily increased to values between 3 and 4, depending on 
tf. Austin at at concluded that the Z model "should only be applied with caution to impact cratering". 
A more general description of crotering flow fields is considered next. 

Implications of a Coupling Parameter 

The cratering process can be conceptually divided into three phases: (1) an early-time phase, where 
the impactor energy and momentum are coupled into the target; (2) an intermediate phase where coupling 
has taken place but the effects of gravity and target strength are small; (3) a late-time phase where 
gravity and strength begin to slow the crater growth. Searches have been made for a measure of "late 
stage equivalence" whereby a single measure of the early- time coupling phase could be used to predict 
the intermediate and late-stage behavior, - Holsopple (1981, 1984) has generalized and quantified this 
idea by introducing a "coupling parameter" which, subsequent to the earlyr-time phase, conveys al! information 
about the impactor., ■ This concept has been used to construct scaling laws for crater ejecta (Ho.usen 
at a/, 1983) and will now be used to investigate the nature of cratering flow fields. 

Suppose that, in the intermediate time regime, the radial velocity of material depends on the impactor 
radius, velocity and density (a, U and 6) and the target density, p, in addition to the position coordinates 
(r,tf) and the time t. The late time regime, where gravity or material strength are important, is not considered 
here. Using a coupling parameter, the radial velocity can be written as 
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r - F( C, p , r, 4; t) 


( 4 ) 


whore P is on unspecified function and the coupling parameter, C, is a scalar-valued function of a, U 
and 6: A dimensional analysis of (4) gives 

f/U « (r/a)‘f *{ d, (r/a)”^“ t Ut/a } (5) 

where 9 is o;i unspecified function and, for convenience, the density ratio' p/fi is assumed constant. 
The exponent fi is related to the dimensions of the coupling parameter and can be predicted from the 
measured dependence of final crater size on initial conditions. In particular, fi should be about 2.5 for 
porous targets (e.g. dry sand) and- 1.8 for nonporous targets (e.g. water, saturated sand, metals). 

In general, the time dependent 7, model as given in eq (1) Is not consistent with a coupling parameter 
because eq (1) implies a product of separate functions of r and t whereas eq (5) does not. However, 
there is a special form of the 7, model which is admitted by a' coupling parameter: the function A(t) 
in eq (1) varies as a power of. time. Consider the three following cases, listed in order of decreasing 
generality, 

1. if A(t) is assumed to vary as a power of t, say t 7 , then the Z model can be viewed as a "local 
approximation" to a coupling parameter. That is, for a small interval of time, the function in eq (5) 
can be approximated locally (in log-log space) by 9 * {(r/a) - * H Ut/a| 7 , 'where y may vary with time. 
Locally, the exponents are related by Z*/J+y(/f-M). Note, if y varies with time then so will Z (fi is 
constant for a fixed material). 

2. If 9 has the special form of u power law, i.e.,‘if y is independent of time, then eq (5) reduces 
to a Z model (again, where A is proportional to t 7 ). 


3. In the special case of steady flow, the Z model corresponds to a coupling parameter with Z»f3, 

Case 1 may explain, at least qualitatively, the observations of Austin et at (1980) that both A and 
Z (in eq 1) increased with time. If, far the case studied by Austin et d, 9 was not a power law and 
the local slope of 9 (i.e., y) started out near zero and increased with time, then Z would start out near 
1.8 ond would increase with time (note, 0~1,8 is expected for their cloy target). A coupling parameter 
may explain the observations of Austin et a/. But without specific knowledge of the time dependence 
observed in their computed flow field it is difficult to perform a rigorous test. 

Piekutowski’s (1980) measurements of particle velocities for a small-scale explosion in dry sand 
are shown below in the form consistent with eq (5). While there is some scatter due to measurement 
error, the data for various "tracer particles" follow a common trend. Hence, a coupling parameter scales 
these data rather well. Piekutowski’s reuslts imply that 9 is a decreasing function of time for dry sand, 
while 9 increased with time In the Austin et a! calculations. This raises the possibility that 9 may differ 
markedly for different materials. 

REFERENCES: Austin, M.G., Thomsen, J.M., Ruhl.S.F., 

0rphaJ,0.L, ond Schultz, P.H.(1 980), Pro c LPS£ 11th. 
23,25-2345. Croft,S.K. (1980), Proe LPSC 11th. 
2347-2378. Grieve.RAF. and Cintafo,M J. (1981), Prec 
i£S£ 22^ 1607-1621. Hdsapple.KA f198l). EOS. 62. 
949. Holeapple,K.A. (19E4) Abstract, this volums. 

Houten,K.R., Schmidt,R.M. and HoltappIs.KA (1983), 
JQR. . 88 f 2485-2499. Maxwell, 0. and Seifert.K.(1974). 
Report DNA 3628F. Maxwell, 0.(1 977) In Impact ami 
Explosion Cr oterino. 1009-1024. Orphal O.L(l977)ln 
Impact ond Explosion Cratering. 907-917. Piekutowski 
AJ.(1980) - Pros IPSC 11th. 2129-2144. 

Schm?dtR.M.( 1983) ■ EOS. 84, 747. SchmWt,R.M.(1984) 

1 0 Abstract, this volume. Schultz, P.H. end Gault, D.E. 

(1979). JGR, 84, 7669-7687. 
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ON CRATER DYNAMICS: COMPARISONS OF RESULTS FOR 
DIFFERENT TARGET AND IMP ACTOR CONDITIONS 

K. A. Hoisapple, University of Washington F5-I0 
Seattle, WA 98195 


A fundamental question of impact cratering mechanics is how the impactor and target 
conditions and material properties affect the crater. This question must be answered in 
order that known results can be extrapolated with confidence to new and different cases. 

Recent studies have uncovered a surprising simplicity in a large number of cratering 
problems. This simplicity is a consequence of a separation of the early from the 
intermediate and late-time mechanics, with a single scalar-valued "coupling parameter" 
depending on impactor radius a, velocity U and mass density d, as C oc aU^ 6 V , being 
the sole measure of the impactor for all of the intermediate to ’ate-time processes. 
Previous applications to final crater geometry (1,2), to ejecta blankets (3) and to crater 
formation time (4) have been reported. 

Here the entire dynamical history of crater growth is considered. Fig. I shows the crater 
depth d versus time t for a wide variety of experimental and calculated craters in non- 
porous targets (5-10), as well as new results for impacts into aluminum. The targets 
include metals, rock, clay and water and the impactors range over a variety of materials 
with the mass densities from 0.01 to 8 gm/cm 3 . The impact velocities range from 2.5 to 
100 km/sec. The results span many decades in size and time. 

• 

Each of these growth curves shows at least some aspects of a generic crater depth versus 
time curve as sketched in Fig. 2, showing three identifiable regimes. In the very early 
time, the flow is one-dimensional and there is no dependence on impactor radius. The 
crater growth proceeds simply as d * vt where v is the initial target particle velocity 
determined by U and the impedances of the two materials. 

As the flow becomes two-dimensional, the impactor radius becomes an important 
parameter. However, after a very short transition time, the sole remaining measure of 
the impactor is the value of the coupling parameter C given above. The exponents and 
V are found to be M = 0.-58 and v = 0.66 for all of the cases shown, arid the crater growth 
can be shown to be governed by the law d<rt M 'U +>1 '. 

Finally, in the late-stages, the growth is inhibited by either material strength Y (in the 
case of small craters) or by gravity g (in the case of large craters; Pgd >> Y ) and the 
scaling of the maximum crater size is governed by the results reported previously. 

These results imply that, in each regime, a different scaling exists that should relate ail 
crater histories. Figures .3, 4, and 5 illustrate this fact showing respectively, the early, 
intermediate and late-time scaling. In each case, the variety of curves do indeed 
superimpose in the applicable regime. 
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K. A. Holajvppie 

1. Hclsappl#,f<A(l981) Coupling parameters In cratering {abstract). In J£5 62. p.949. 

2. Holsopp!e,K.A.(l9S3) On the existence and Implications of coupling parameters In cratering mechanics 
(abstract). In Lunar and Planetary Science XIV.sp.3 19-320. 

3. Heui«n.K.R.^chmidh9.U..Holsappl>.KJi.( 1 982)J. Geoohvs. Res. S3t'B3) . pp.2485-2499. 

4. Schmfdt,R.M.(1981) Sealing crater time of formation (opsrrcct). In COS 62.No45 .aS44. 

5. Dienes, J.K.,Wdsh,J.M.(1 970) In High Veleclfv Imocef Phenomena ,, ed. by R.KInslow, Academic Press, New 
York, pp45-104. 

6. Gault.D.L.Scn net,C.P.( 1 95 1 ) In Geological Implications of Impacts of Large Asteroids and Comets on the Earth. 
ed.by LT.Si!ver arid P.H.Schultz. 

7. Klneke,J.H.,(196l) Observations of crater formation In ductile matsrlair. In Proceedings of the Sth Svmcoslum 

5. O'Keefe, J.D., Ahrens, TJ.(1 982) Cometary and 

Meteorite Swarm Impacts on Planetary Surfaces. J. 
Geophvs. Res. 

9. Crphcl.D.L, Borden .W.F., Lcrson.SA, 

Shultz, P.H.(1980) Impact Melt Generation and 
Transport. In Proc. tlth Lunar end Plcnetory 
Science Conf.. pp.23C9-2323. 

10. Austin, M.G., Thomsen, J.M., Ruhl,S.F., 0rphol,0.L, 
Shultz, P.H.{ 1980) Calcuiattonal Investigation of Impact 
cratering dynamics: material motions during the crcrer 
growth period. In Pr oc, 1 1 th Lunar end Planetary 
Science Conf. .pp.2325-2345. 


On Hvoervelocitv Imoact. voU. part 2, pp339-370. 









tr 

s 






ORIGINAL PAGE »« 
OF POOR QUALITY 


FEASIBILITY OF DETERMINING IMPACT CONDITIONS FROM TOTAL CRATER MELT 

Michael D. Bjorkman, The Boeing Company, Seattle, VA 98 124 

In the past, authors have used parameters such as kinetic energy or momentum to scale 
crater volume and totai melt volume (sum of that in crater and ejecta) produced by 
impact. The converse procedure; determining impact characteristics from, for example, 
the crater volume, will not determine impactor radius, a, and velocity, U, individually, but 
only as a prodict of powers of a and U. Recently, Grieve and Cintaia (1) proposed to 
solve crater volume and total meit volume relations simultaneously to determine a and U 
individually for a terrestial crater. Their analysis pointed out that published total melt 
volume scaling relations varied from energy scaling (2) to momentum scaling (3), with no 
indication of the domain of applicability. The coupling parameter notion, described 
below, has the potential to unify the melt scaling relations and identify which materials 
nearly energy scale and which momentum scale. 

The author examined the scaling of depth of melted planetary material along the axis of 
impact and showed (4) that a coupling parameter of the form C * a U M d v , u* 0.38, V * 
0.66, which governs crater depth in nonporous materials like rocks and metals, also 
governs melt depth. The parameter 6 Is impactor density. This result necessarily follows 
from the fact that C determines sources which produce equivalent flow fields at 
intermediate and late times (3). 

If the source variables a,,U and 6 only enter the scaling of melt depth D m via C, then, 
D m * f(C, p , E m ), where p is the planet density and E m is the specific energy of the 
Hugoniot state from which isentropic release ends at the 1 atm point on the liquidus. 
Dimensional analysis results in D m /'a * (U 2 /E m ) u/ * C6yp) v . Figure 1 is a plot of meit 
depths calculated with hydrocodes (2,6,7). If the above scaling relation holds, then the 
curve should be horizontal. Figure 1 shows that C does govern the melt depth for 
v 2 /E m - 8 (where v is the particle velocity at impact). 

A similar analysis was performed for total melt volume (2,6-8) but did not demonstrate 
the coupling parameter governs melt volume. Figure 2 shows the 'data are not horizontal, 
and appear to be closer to energy scaling, p = 0.67. Perhaps crater and melt volume 
scale differently, however, this is difficult to understand, because they both stem from 
the same flow field, which has been shown to scale with Us 0.38 (9). The different 
scaling for the melt depth and melt volume is being looked into, 2 nd is suspected to result 
from free-surface effects. 

The meit depth *nd melt volume scaling relations are both significantly different from 
momentum scaling, p = 0.33, however. The momentum scaling result (3) for total meit 
volume was obtained from an analytic model which prohibited ejection of material from 
the crater. This prevents multiplication of the downwards directed momentum and thus 
precludes any result other than momentum scaling. Momentum scaling of crater volume 
has been observed only for porous materials; it is therefore possible that the scaling 
relation from ref (3) applies to porous materials only. 

If the present approach is correct, then it is not possible to distinguish the effects of the 
impactor radius and its velocity U separately using late time phenomena like totai crater 
melt, crater dimensions or phase transitions in rocks below the crater. Only material 
from near the impact would retain information about a and U separately. Search for this 
material is complicated by ejection to great distances. Even if crater meit volume does 
scale with energy, U = 0 .67 and crater volume with M = 0.38, the difference is not great 
enough to allow one to determine a and U separately to any useful accuracy. 


FEASIBILITY OF DETERMINING IMPACT CONDITIONS FROM TOTAL CRATER MELT 

8 ►— *— #***> * 

Michael D. Bjorkman 


t. Grieve, R, A. F. and Cintala, M. 3. (1981) Proc Lunar Planet Sc: 12th . p. 1607-1621. 

2. O'Keefe, 3. D. and Ahrens, T. 3. (1977) Proc lunar Scl Con* &h . 573337-3374. 

3. Gault, D. E. and Heitowit, E. D. (196i) in Proc 6th Hypcrvelocitv Impact Svmp. 2 
p. 419-456. 

4. Bjorkman, M. D. (1983) EOS (abstract), 64 (43) p. 747. 

3. Bjorkman, M. D. and Hoisapple, K. A. (T383). Proceedings of 3rd APS Conference on 
Shock Waves in Condense Media . 

6. Hoisappie, K A, (1984) On crater dynarmi'4: Comparisons of results for different 
target and im pact or conditions. This volume. 

7. Orphal, D. L. et ai (1980) Proc Lunar Planet Sd 1 1th. p. 2309-2323. 

8. Bjork, R. L. et ai. (1967) f?inaiytical Study oi Impact Effects as Applied to Meteoroid 
Hazard . NASA report CR-/#. 

9. bienes, 3. K. and Walsh, 3. M., in High Veiocity Impact Phenomena ed. by R. Kinsiow 
(Academic Press, New York, 1970) p. 46-104. 
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Figure 1. Scaled meit depths. 



Figure 2. Scaled melt volumes. 
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